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We present the POLEMIC A^ (Polarized Expectation-Maximization Independent Component 
Analysis) algorithm which is an extension to polarization of the SMICA^ temperature com- 
^ ■ ponent separation method. This algorithm allows us to estimate blindly in harmonic space 

\^ ' multiple physical components from multi-detectors polarized sky maps. Assuming a linear 

, noisy mixture of components we are able to reconstruct jointly the electromagnetic spectra 

(f) ■ of the components for each mode T, E and B, as well as the temperature and polarization 

spatial power spectra, TT, EE, BB, TE, TB and EB for each of the physical components 
r*"** ■ and for the noise on each of the detectors. This has been tested using full sky simulations 

of the Planck satellite polarized channels for a 14-months nominal mission assuming a simple 
pH . linear sky model including CMB, and optionally Galactic synchrotron and dust emissions. 

Or 

2 ' 1 Introduction 

Mapping the Cosmic Microwave Background (CMB) polarization is one of the major challenges 
of future missions in observational cosmology. CMB polarization is linear and therefore can 
be described by the first three Stokes parameters I, Q and U which are generally combined 
to produce three fields (modes), T, E and B. The polarization of the CMB photons carries 
extra physical informations that are not accessible by the study of the temperature anisotropies. 
Therefore its measurement helps breaking down the degeneracies on cosmological parameters 
as encounter with temperature anisotropies measurements only. Furthermore, the study of the 
CMB polarization is also a fundamental tool to estimate the energy scale of inflation. 

However, CMB polarization is several orders of magnitude weaker than the temperature 
signal and therefore, its detection needs an efficient separation between the CMB and the astro- 
physical foregrounds which are expected to be significantly polarized. 

A direct subtraction of these foreground contributions on the CMB data will require an 
accurate knowledge of their spatial distributions and of their electromagnetic spectra. But these 
latter are not yet well characterized in polarization. 

To try to overcome the above limitations, a great amount of work has been dedicated to 
design and implement algorithms for component separation which can discriminate between 
CMB and foregrounds. We present here the PolEMICA (Polarized Expectation-Maximization 
Independent Component Analysis) algorithm which is an extension of the Spectral Matching 
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Independent Component Analysis (SMICA) which has been developed to consider both a fully 
blind analysis for which no prior is assumed and a semi-blind analysis incorporating previous 
physical knowledge on the astrophysical components. This extension allows to estimate jointly 



the temperature and polarization parameters from a set of multi-frequencies I, Q and U sky 
maps. 



2 Model of the microwave and sub-mm sky 



To perform the separation between CMB and the astrophysical foregrounds, the diversity of 
the electromagnetic spectra and of the spatial spectra of the different components is generally 
used. Observations from a multi-band instrument, for the Stokes parameters I, Q and U, can 
be modeled as a linear combination of multiple physical components leading to what is called a 
Multi-Detectors Multi-Components (MD-MC) modeling. 

Assuming an experiment with n u detector-bands at frequencies Vi and n c physical compo- 
nents in the data, working in the spherical harmonics space, we can model the observed sky for 
X = {T,E,B}, for each frequency band and for each {£, m} 
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where y^ m is a vector of size (3- n v -n£ ■ n m ) containing the observed data, s c £ ' m is a (3 • n c • ri£ ■ n m ) 
vector describing each component template and is a vector of the same size than 
accounting for the noise. A" ,x is the mixing matrix containing the electromagnetic behaviour 
of each component and is of size (3 ■ n u ) x (3 • n c ) 

The aim of the component separation algorithm presented in here is to extract A v ^ , s^ m 
and n£ m from the y^ sky observations. 



3 A MD-MC component separation method for polarization 

To reduce the number of unknown parameters in the model described by equation (pQ), it is 
interesting to rewrite this equation in terms of the temperature and polarization auto and cross 
power spectra and to bin them over £ ranges. 



Ry(b) = AR s {b)A A +Rn(b) (2) 

where R y (b) and R n (b) are (n u • 3) x (n v ■ 3) matrices and R s (b) is a (n c • 3) x (n c • 3) matrix. We 
assume that the physical components in the data are statistically independent and uncorrelated 
and that the noise is uncorrelated between channels. 

To estimate the above parameters from the data we have extended to the case of polarized 
data the spectral matching algorithm developed in smicaS for temperature only. The key issue 
of this method is to estimate these parameters, or some of them (for a semi-blind analysis), 
by finding the best match between the model density matrix, R y (b), computed for the set 
of estimated parameters and the data density matrix R y {b) obtained from the multi-channel 
data. The likelihood function is a reasonnable measure of this mismatch. We have extended 
this method to jointly deal with the temperature and polarization power spectra and also to 
estimate the TE, TB and EB cross power spectra^. 

The maximization of the likelihood function is achieved via the Expectation-Maximization 
algorithm (EM)^ . This algorithm will process iteratively from an initial value of the parameters 
following a sequence of parameter updates called 'EM steps'. By construction each EM step 
improves the spectral fit by maximizing the likelihood. For a more detailed review of the spectral 
matching EM algorithm used here, see^. 



4 Simulated microwave and sub-mm sky as seen by Planck 



Following the MD-MC model discussed above and given an observational setup, we construct, 
using the HEALPix pixelization scheme^ and in CMB temperature units, fake I, Q and U maps 
of the sky at each of the instrumental frequency bands. For these maps we consider three main 
physical components in the sky emission: CMB, thermal dust and synchrotron. Instrumental 
noise is modeled as white noise. 

The CMB component map is randomly generated from the polarized CMB angular power 
spectra for a set of given cosmological parameters. In the following we have used Hq = 71 km • 
s- 1 ■ Mpc~\ n b = 0.044, n m = 0.27, ^ A = 0.73 and r = 0.17 that are the values of the 
cosmological concordance model according to the WMAP 1 year results^. 

For the diffuse Galactic synchrotron emission we use the template maps in temperature and 
in polarization provided by^. Here we have chosen to use a constant spectral index equal to the 
mean of the spectral index map, a = —2.77, so that the simple linear model of the data holds. 

In the case of the thermal dust we dispose of few observational data of the polarized diffuse 
emission and to date no template for this is available. Thus, we have considered a power-law 
model, renormalized to mimic at large angular scales the TE cross power spectrum measured 
by Archeops at 353 GHzl i", Q and U full-sky maps are generated randomly from these power 
spectra. We extrapolate them to each of the frequency of interest by assuming a grey body with 
an emissivity of 2. 

Noise maps for each channel are generated from white noise realizations normalized to the 
nominal level of instrumental noise for that channel. 

We have performed sets of simulations of the expected Planck satellite data to intensively 
test the algorithm presented above. We present here results from 300 realizations considering 
full-sky maps at the LFI and HFI polarized channels, 30, 40 and 70 GHz for LFI and 100, 143, 
217 and 353 GHz for HFI for a nominal 14-month survey. We have simulated maps at n s jd e 
= 512 which permits the reconstruction of the angular power spectra up to i ~ 1500. The 
reconstructed spectra will be averaged over bins of size 20 in t. 

5 Results 

We have applied the PolEMICA component separation algorithm to the simulations presented 
above. From them, we have computed the data density matrix R y and applied the algorithm. 
We simultaneously estimate the R s , R n and A matrices, with no priors, for temperature and 
polarization. To ensure the reliability of the results we have performed 10000 EM iterations and 
checked, for each simulation, the convergence of the EM algorithm. 

We present in figure [U the reconstructed CMB power spectra. We can see that for the TT, 
TE, TB and EB spectra, we are able to reconstruct the Cg over the full range of i values that are 
accessible at this pixelization resolution (^ max ~ 1500). The EE spectrum is recovered accurately 
up to I ~ 1200. For smaller angular scales, a bias appears. This bias is a pixelization problem 
that would occur at a larger t if the resolution was higher. The BB spectrum is reconstructed 
up to i ~ 70. For larger £, the reconstructed spectrum is residual noise arising from the fact 
that the convergence of the EM algorithm is slow and therefore we have not properly converged. 
This bias appears in our separation when the signal over noise ratio is below 10~ 2 and does not 
affect the reconstruction of the other parameters. Even if we were able to avoid this effect, the 
recovered BB spectrum would be compatible with zero for i > 70 thanks to the size of the error 
bars. 

The power spectra from our input synchrotron and dust emissions are recovered with effi- 
ciency up to I ~ 1500 for TT, EE, BB, TE, TB and EB^. Power spectra of the noise in 
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Figure 1: Reconstruction of the CMB power spectra for C/ 1 , Cf , Cf , Q , and Cf" at 100 GHz in 

mKcmb performed on Planck full sky maps simulations. Crosses represent the reconstructed spectra, solid lines 
the input model. Error bars are the dispersion over ~ 300 simulations. 



temperature and in polarization are also fully reconstructed^. 

The mixing matrix A elements corresponding to CMB and dust emission are recovered 
efficiently, for temperature and polarization. For the synchrotron emission, mixing matrix ele- 
ments corresponding to polarization are well recovered and those corresponding to temperature 
are biased at intermediate frequency values^. This bias is due tu a slight mixing up between 
synchrotron and CMB in temperature. It does not happen in polarization where the synchrotron 
dominates the CMB. This bias can be avoided by the adjunction of priors in the separation, like 
for example assuming an equal electromagnetic spectrum in temperature and polarization for 
each component^ 

To evaluate the impact of foregrounds in the determination of the CMB temperature and 
polarization power spectra we have compared the results of the presented analysis to those on 
simulations that contain only CMB and noise. In the presence of foregrounds, the error bars on 
the reconstruction of the CMB power spectra are increased by at least a factor of two both in 
temperature and in polarization^. Therefore, although the foreground contribution in the data 
can be removed, it significantly reduces the precision to which the CMB polarization signal can 
be extracted from the data. 
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